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Abstract 

The Self-Consistent RPA (SCRPA) equations in the particle-particle 
channel are solved without any approximation for the picket fence 
model. The results are in excellent agreement with the exact solutions 
found with the Richardson method. Particularly interesting features 
are that screening corrections reverse the sign of the interaction and 
that SCRPA yields the exact energies in the case of two levels with 
two particles. 
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1 Introduction 



The approach to the many body problem is often divided into the evaluation 
of a hierarchy of correlation functions. First comes the one body or mean 
field theory. Second the two body correlation functions, and so on. There 
is a large consensus on the mean field or Hartree-Fock level so that practi- 
cally in any field of theoretical physics one basically understands the same 
by this approach in spite of the fact that there may be differences in detail. 
However, already on the next level of complication, i.e. considering two body 
correlations, the strategies in various fields start to diverge strongly. Let 
us mention first a few of them: Brueckner-Hartree-Fock p|, Jastrow ansatz 
for correlated wave functions, the hypernetted chain method, correlated ba- 
sis functions 0, coupled cluster theory |3J, Gutzwiller ansatz for correlated 
ground states [§J, and quite a few more are examples of different strategies 
to tackle the problem of strongly correlated quantum systems. There seems 
no common consensus about which of the strategies is best, all seem to be 
more or less tailored for specific problems. In this situation it seems legiti- 
mate to explore further routes which have been not sufficiently exploited in 
the past. It is the latter case which we are investigating since a couple of 
years, further developing the so-called Self- Consistent RPA (SCRPA) the- 
ory which has been initiated mostly in nuclear physics with the pioneering 
work of Hara ||, and further developed by Rowe ||. Independently and 
using different techniques the same approach has been developed in plasma 
physics where it is called cluster Hartree-Fock theory [f?j]. Very early quite 
similar ideas have also been put forward by the Japanese school of Marumori 
et al P| . The recent theoretical developments along this line and the various 
successes in applications to model cases [[], |IJj, [IT], [12], [IB], [IJ, |IJ| |IB[ make 
us firmly believe that there exists a great potentiality in this approach which 
is very general and can be easily adapted to practically any kind of many 
body problem where two body correlations or even higher correlations are of 
importance. 

The SCRPA scheme has indeed, on the level of the two body correlations 
or fluctuations characterized by a correlated pair of Fermion operators like 
Ql — J2kk' ^kk' a k a k' or Sfcfc' ^kk ,a k a k', aspects quite analogous to ordinary 
Hartree-Fock theory characterized by the single operator = X^^fc* 3 ! ■ As 
we will show in the main text, standard Hartree-Fock theory and SCRPA 
or Cluster Hartree-Fock theory can be derived from a common variational 
principle leading to the equations 
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{6q, [H,ql]}) = e a ({5q,ql} 



8Q,\H,Qt]\) = E v (\8Q,Qt\) , (2) 



where {, } is the anticommutator and [, ] the commutator, and (• • •) stands 
for the mean value with respect to vacua defined by q a |0) = or Q v |0) = 0. 
Equations (|]) and (0) are not only formally very similar. Also their physical 
content is quite analogous. Whereas (jl]) describes the motion of one particle 
in the average field created by all the others, the second equation (g) describes 
the motion of a quantal fluctuation (a correlated Fermion pair) in the average 
field created by all the other fluctuations (pairs). A nucleus, for example, 
can be considered in some approximation as a gas of independent zero point 
vibrations (elementary excitations) as well as a gas of independent particles. 
Both gases create their own mean field and in general they will be coupled. 

As with any theory for correlation functions also SCRPA is very diffi- 
cult and demanding with respect to numerical implementation in the general 
case. However the recent solution of the many level pairing model is suffi- 
ciently complex to state that with respect to a brute force diagonalization, 
the SCRPA gives a tremendous simplification of the solution of the many 
body problem with very good accuracy which, at least for the case consid- 
ered, does not fail in any qualitative respects. It is our intention in this work 
to take up again the SCRPA solution of the many level pairing model of Ref. 



15| in more detail. In first place we circumvent an approximation used in 
15f so that in this work the SCRPA scheme is carried through in full without 
any approximation. On the basis of this achievement we will then be able to 
discuss the advantages and eventual shortcomings of the method. We also 



will be able to assess the deficiency of the approximation we adopted in |T5 



as well as other more severe simplifying reductions of the full SCRPA scheme 
such as the renormalized RPA (r-RPA) which, because of its relatively simple 
numerical structure, has become quite popular in the recent past [17|, [R| [L9| . 

In detail the paper is organized as follows: in Section 2 the picket-fence 
model is introduced, in Section 3 it is reformulated to make explicit its 
particle-hole symmetry, in Section 4 the SCRPA formalism is presented, in 
Section 5 the occupation numbers and correlation functions needed to close 
the equations are constructed, in Section 6 the approximations used to obtain 
simpler RPA formalisms are discussed, in Section 7 sum rules related with the 
conservation of the particle number are presented. The explicit construction 
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and the problems associated with the SCRPA ground state are discussed in 
Section 8, numerical results are shown in Section 9, and the conclusions in 
Section 10. Some useful mathematical results are given in the Appendix. 



2 The Picket-Fence Model 

The picket fence or multilevel pairing model has been introduced by Richard- 
son in 1966 [Q to describe deformed nuclei. The model has the great advan- 
tage of being non-trivial (in the sense that it can not straightforwardly be 
diagonalised for an arbitrary number of levels) and still being exactly solv- 
able (by some kind of Bethe ansatz). However, in nuclear physics, aside from 
some rare considerations pi] , the model has not been exploited very much 
probably because it has been judged too crude for the description of real 
nuclei. Still the model contains very interesting physics and it has recently 
been revived in the context of ultra-small superconducting metallic grains 



22 . One of the most interesting aspects of the model is that the exact solu- 
tion reveals a transition between the superfluid (or superconducting) regime 
and the normal state which is completely smooth, i.e. no sign of any abrupt 
phase transition from one state to the other can be detected as a function of 
the system parameters \ T5 \. 



The original Hamiltonian of the model is given by 

H = J2(e l -\)N l -Gf:P?P 3 , (3) 

i=l i , j = 1 

with 

N r = c\ Cl + cl t c. h P} = c\cU, (4) 

where cj creates a particle in the i-th level with j — \ and m = | and cLj 
with m = — |. Q is the total number of levels, G is the pairing interaction 
strength and £j = % e. The chemical potential A will be defined such that the 
Hamiltonian preserves particle-hole symmetry. Each level has degeneracy 
two and we will assume that the system is half filled with number of pairs 

n = n/2 . 

The particle (p) and hole (h) states are defined by 

N h \HF) = 2, N P \HF) = , (5) 
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where \HF) is the ground state of the Hamiltonian (|3|) with G=0. The 
particle states p correspond to e p > A and the holes h to e p < A. 

In this case, with no partial occupations allowed, the following relation is 
fulfilled: 



PjP l + P l P} = l, 



which implies 



N> = 2P}Pi 



(6) 
(7) 



3 The particle-hole symmetry 

The commutation relations between the operators defined in (2) are 



p. pi 



Sij (1 - Ni) , 



N t ,P] 



-28ijPj 



To make explicit the particle-hole symmetry we will make the following 
particle- hole conjugation 



Cp dp , Cp dp , 



di 



ck = -4- 



The new operators M, Q and in terms of d' (d) are 



2 - M h , N n = M n 



Hi — * — lvl hi — 

Their commutation relations are 



r i, 



-Qt. 



Pi = Ql ■ 



(9) 



(10) 



Qh, Qh' 



5 pp/ (l-M p ) , 
= 8 hh , (1-M h ) 



M p ,Ql] =25 P p,Ql 



v ' 

M h , Ql] = 25 hh ,Q ] 



h ■ 



[M p , Q p i] — —25p p iQp , 
[M h ,Q h ,] = -25 hh ,Q h .(11) 



The single particle energies are e p = e(N + p) and = e(N — h + 
1), with p, h = 1,...,N , and = Particles (p) and holes (h) are 

numbered starting from the level closest to the Fermi level. We use the 
chemical potential to restore the particle-hole symmetry: 



A = e { N + 



G 
~2 



;i2) 
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With this definition the Hamiltonian reduces to 



H = -eN 1 



N 

E 

p=h=l 



G 

e \P-~ ) +77 



{M p + M h ) 



G E QlQv - G E QlQ/>' + c E (QjQi + QpQh) ■ (13) 



yjy/ 



In this form the complete symmetry between particle and hole states 
becomes evident. This will greatly facilitate later the formal and numerical 
aspects of our theory. 



4 The RPA formalism 

The basic ingredients of the SCRPA approach in the particle-particle channel 
are the two particle addition operator 

4 = E^^-E^Q,, (14) 

p h 

and the removal operator 

R\ = -Y, Y pQp + Y, x tQ{, (is) 

P h 

where Q p = Q p /y/l - (M p ) and Q h = Q h /y/l-(M h ). 

Following Baranger in his derivation for the single particle mean field 
equations J23J, we define the following mean excitation energy 



fy. = {E a{N+2) ( E N + 2 ~ E° N+2 )\(a\Al\0)\ 2 

+ E/3(iV-2) {En-2 ~ E N~2) I (01 A,|0)| 2 

+2/iW EJH4|0}| 2 - 2 M W Ey3 |(/3K|0}| 2 ]} / (16) 

{E Q(J v +2) K«I4|0)| 2 - E^- 2 ) |<M,|o)| 2 } , 

where 2/iW = (±)l/2(£'^ r+ / 2 — E 1 ^-) are the chemical potentials and 

are the in principle exact eigen values of H of Eq.(l) and |a), |/3), |0) are the 

corresponding exact eigenstates. 



7 



Expression ( ]TE| ) can be considered as an average excitation energy taking 
into account the spectra of both the N + 2 and N — 2 systems. In fact 
represents the energy weighted sum rule over both spectra divided by the 
non-energy weighted sum rule. 

Indeed expression (|16D can also be written as 

(0P„[ff,4]]|0) 

A* ~~ /r.1 r A aU lr>\ ' \ ' 



which makes the analogy with the usual energy weighted sum rule in the 
particle-hole channel quite obvious. The terms in ( |i6| ) involving the chemical 
potentials are needed in order to give the correct origin of the energy 
spectra. This is best seen in supposing that = which leads to 

n _ 2 _ EaQEg+g - n +2 )|(«l4|0)| 2 + Ep(E p N _ 2 - n_ 2 )|(/3|^|0)| 2 
" ^ (0|[A„4]|0> 

(18) 

It means that the energy weighting starts from two times the chemical po- 
tential, as it should be. 

The first step is now to minimize the average two particle energy ( |i6|) 



with respect to the amplitudes X, Y . This leads straightforwardly to the 
following set of equations: 





(19) 



where 



A pp/ = (0\[Q p , [H,Ql,]}\0) 

V \ 2 (e(p - i) + § ) + ^^-yUE Ql - E Q*)Q,> 

_ G _ ((1-M P )(1-M P 0> 



^(1-(M P ))(1-(M P ,)) 



Ju=<oip? > ,[g,gt]]|o> = G ff"^"^, . ( 2 °) 

V /(1-(M P ))(1-(M,)) 



cw = (o|[Q„[#,Ql]]|o> 
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>hh> 



-^<h-l) + ^-2^ ){ Ql { -^Qt i + ^Q hl)) 

{{I - M h ){\ - M h ,)) 



+G 



y/(l-(M h ))(l-(M h ,)) ' 



Due to the particle-hole symmetry, the removal mode satisfy exactly the 
same equations, implying that both modes have exactly the same excitation 
energies and wave functions. To advance this conclusion we are assuming the 
following relations: 

(M p ) = (M h=p ) , 

(QlQv) = (Ql= P Q h >= P '), (QhQp) = (QlQl) (21) 

(M p M p/ ) = (M h=p M h , =p ,), (M h M p ) = (M h , =p M p , =h ), 

which are shown below to be consistent with the above equations. 
It means that 

Xg = ±X^, YC = ±F p A = T. (22) 

The forward and backward amplitudes X, Y fulfill the normalization condi- 
tions 

V h 

E^'-EW = ^v, (23) 

h P 

P h 

and the closure relations 

A* A 

]TX^-£y^ = <W , (24) 

E^ P A -EW = o. 

A At 

The expectation values of the commutators are 
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R\, Ry 



8, 
5„ 



4t Z?t 



W > 

([A fl ,Ry})=0 , 

Rx,A\ 



(25) 



. (26) 
With the help of these equations we can now invert fll4 ,15[) which yields 



Q 



(M p ) 



Q h = y/l- (M h ) 



E^4 + E^a 

, V A 

E*^+En M 4 



L A 



(27) 



An important point to recognize at this step is the fact that the above 
SCRPA equations can also be derived using the equation of motion method 
advocated for example by Rowe . A basic ingredient to this latter method 
is that one assumes the existence of a vacuum state |0) = \SCRPA) such 
that 

A^ \SCRPA) = R x \SCRPA) = 0. (28) 

Therefore Eq. (p5[) must be considered as an inherent additional relation 
belonging to the whole SCRPA scheme. The expectation values in ( ^0|) shall 
then be evaluated with this ground state and, as we will show in a moment, 
this will permit to entirely close the system of equations. 



5 Closing the system of equations 

The aim is now to express the expectation values which figure in the RPA- 
matrix (P0| ) entirely and without approximation by the amplitudes X, Y. 
Using ([27], 28|) this can easily be achieved for most of the expectation values. 
Indeed one directly verifies: 



(QlQ P >) 



<M P ))(1-(M P ,))EW 
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(QhQh') = J(i ~ {M h ))(i ~ (%))EW. (29) 



(Q h Q P ) = (QlQl) = V(i - (M h ))(i - {M p ))Y,Y p x xl 

A 



Using (HJIOD it is direct to show that 

(M„) = 1 — — . (M h ) = 1 — — , (30) 

v p/ i + 2E(>; A ) 2 1 + 2EW) 2 ' 1 ; 

which together with ( |22| ) implies 

(MpJ) = (M h=l ) , (31) 

reflecting again the particle-hole symmetry. Using (|22|) and ([31]) it is also 
direct to show that 

(QlQ P >) = (QU p Q h ^). (32) 

These are the first three of the five equalities advanced in (|2T1) . 

Knowing these expectation values we can evaluate the SCRPA ground 
state energy: 



N 



(H) = -eN 2 + £ 



p=h=l 



G 

£{P- o ) + 77 



«M P ) + <M h » 



-GY,{QlQ*) - GJ2(QlQh') + ((QlQl + QpQhi) • (33) 

pp' hh' ph 

The SCRPA correlation energy is 

E^ RPA = (H)+eN>. (34) 
For comparison, the RPA correlation energy is 

eZ a = -Y.e»Y.K\ 2 - (35) 

n p 

In order to fully close the set of SCRPA equations we still must express 
the correlation functions (MiMj) through the RPA amplitudes. This can 
also be done exactly, though it is somewhat involved. For this reason we 
approximated (MiMj) by (Mi)(Mj) in In what follows we now present 
the full derivation. 
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Using (|j) and (|7]) it is direct to demonstrate that 

NiNi = 2N U (36) 

which implies that 

M p M p = 2M P , M h M h = 2M h . (37) 
It is also simple to demonstrate that 

N t N 3 =AP}p}P 3 P i for i^j, (38) 
which implies the three sets of equations 

M p M p ,= AQlQl,Q p ,Q p ioip^p', 
M p M h = M p + M h -2 QlQhQlQp - 2 Q{QpQ\Qh , (39) 
M h M h t = ^QuQh'Qh'Qh for h^h' . 

We have then Q expectation values which are known: 

(M p M p ) = 2(M p ), (M h M h ) = 2(M h ). (40) 
For p 7^ p',h 7^ h! the equations are 

(M p M p/ ) = 4(1 - (M p ))(l - (M p/ )) E E Y p ^'Y p )^(R x R Xl R{A') - 

AA'AiA 2 

(M p M h ) = (M p ) + (M h ) (41) 

fpY p X^ 1 X^ 2 (R x R Xl R\ 2 , 



-2(l-(M p »(l-(il4»E £ Y p %*X^X^(R x R Xl R{R{i) 



XX' AiA 2 

-2(1 - (m p »(i - (M fc » e e wiri^v^^i') > 

(M h M h /) = 4(1 - (M h »(i - (M h /)) e E W'y«y;(v w 44) • 

MM' WW 

Using the identities given in the Appendix we build up the system of 
equations 

£ <M *' Mfe)(sS ' " 4(i- ( Awi*a 2 -(A4 2 )) 2) + 

Y.(Mk 1 M pi )C^ 1 + (42) 
y E{M K M pa )D** a = F hh ,, 

PlP2 
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£<M fcl M ha )i# ha + 

g(M fa M p ,)(C & , + 2(1 _ (M ^_ (Mpi))2 ) + (43, 
^2(M Pl M P2 }DpP p2 = F hp , 

E<^i^pi>C^+ (44) 

V7M fl/f \(T)PP' ^ppi Vpg \_ F 

^{M^M^U^ 4(1 _ (Mpi>)2(1 _ (Mp2>)2 J 

were we have introduced the auxiliary matrices 

^ = -2E P1 (1 - (M Pl ))y 2 pP X' Pl + 2E fcl (l - (M fcl »^^ fcl 

-E fclfca 5«f fca (l-<M fcl >-<M fca » (45) 

-e, 1p1 ^; 1 (i-(m p1 )-(m, 1 )) 
-e p1P2 ^; 2 (i-(m p1 )-(m P2 )) , 



CZ' P1 = -(y^+^PiW. (46) 

Dpip 2 = (3^ PP iD^ P ' P2 ) + 3^ PP i3V P2 3V P i3^ PP2 • 

% = «M fc > + (M p »/(2(1 - (M,)) 2 (l - <M P » 2 ) 
— 2E Pl (l — (M Pl )) (y ppi Z pih — X ppi Zl pi 

-2E^(i - (-M/u)) (^^y^j - z^ hi x^ hi 

-E hlh2 Bf ih2 (l-(M hl )-(M h2 )) (47) 
-E, 1P1 < P1 (1-(M P1 )-(M, 1 )) 
-E PlP2 ^f P2 (l-(M Pl )-(M P2 )) , 
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Bhih 2 ~ {yh^Zh-ip) 2 + yhh 1 Zh 2 pZh lp yhh 2 

+ (Xhh 1 Zph 2 ) 2 + Xhh 1 Zph 2 Z p h 1 Xhh 2 i 

Cfefpi — -(y ppi Xhhi + Z Plh Z phl ) 2 (48) 

(3^/i/ii^ppi 4" Zf ll p2f ipi ) , 

~^pip 2 ~ (z P1 hy PP2 ) ~i~ z pl i l ypp 1 Zp 2 i l y P p 2 

J ^i.Z} l p 1 Xpp 2 ) + Z} l p 1 Xpp 1 Zf ip2 Xpp 2 . 



F w = 2E P1 (1 - (M pi ))2% pi Zh, pi - 2£ fcl (l " (Mk,))^^ 

- E/nha BjJUl " (M fcl > - (M, 2 » (49) 
-E ftipi Cg 1 (l-(M pi )-(M, 1 )) 
-E P1P2 ^ 2 (1-(M P1 )-(M P2 )) , 



Cftfpi = -(-^ft'pi^/u + z hpi y h , hl ) 2 , 

Dpip 2 = (Zh'p 1 Z} l p 2 ) 2 + .2^/ p l Zf ip2 Z] ipi Zf l i P2 ; 



with 



^feftl — r — — — , <^-ppi — l^i p 



V v /(l-(M,))(l-(M hl ))' m ^Va-WXi-^))' 



A 



v /(i-(M h ))(i-(M p ))' P W^^wd-W)' 



Due to the symmetry (M h M p ) = (M p M h ) there is no need to include an 
explicit equation for (M p M h ). 
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Notice that the above matrices have the following symmetry properties 
nw' — tdp'p nPP 1 — np'p r>PP' — Dp'p 

D h 1 h 2 ~ n h 2 h 1 i Wipi — °/upi> U pip 2 ~ U P1PV \ 0Z 1 
r>hh' ryh'h /^thh 1 r^th'h, -pMi 1 r^h'h 

n hih 2 - n h 2 hv U h lPl - ^hipi) U VW1 ~ JJ P2Pi- \ 06 ) 

They imply 

Fpp> — F p ' p , F hh i = F h i h , (54) 

and guaranty that (M p M p >) = (M p ,M p ) and {M h M h ,) = {M h ,M h ). 

Due to the particle- hole symmetry, Eqs. ( p2|) and (|3T|), these matrices 
possess also the following symmetries 

Xhh' = Xp=hp'=h', ypp' = yh=ph'=p', Z p h = Zh'=pp'=h- (55) 

which reflects in the A,B,C,D matrices as 

r>hh' T-ip=hp'=h' j-^hh' f>P=hp'=h' 

D h\h 2 ~ p 1= hip 2 =h 2 i ■ U pip 2 ~ h 1 =pih 2 =p 2 i 

C hi P1 = C h7pi P ~ k ' F PP' = F h=ph>=p'i (56) 

thus implying that the particle-hole symmetry is fulfilled in the following 

sense: 

(M P M P ,) = (M h=p M h , =p ,), (M p M h ) = (M p , =h M h , =p ). (57) 

These are the last two equalities which were assumed valid in (EHf). They 
show that the whole formalism is consistent with the particle - hole sym- 



metry. With these symmetry relations Eqs. (|42||4^J44]) together with Eqs. 
(^5],^B],^7|,^5|,PD|,|5T1D can be solved numerically. 



6 Simpler RPA formalisms 

In a previous work the SCRPA equations (^Up were solved under the approx- 
imation 

{MMj)- {Mi){Mj). (58) 

This replacement close the SCRPA equations, and will be denoted SCRPAi 
formalism from here on. 

The renormalized RPA (r-RPA) formalism allows the consideration of 
the Pauli principle in a simpler way. The commutation relations between the 
operator Qi, Qj are replaced by their expectation value, i.e. 

Qi, Qj] -«y fi (l-<M i ». (59) 
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The RPA equations take the much simpler form 

4? PA = 2<V (c(p - |) + f ) - G^(1-(M P ))(1-(M P ,)) 
^, RPA = Gy/(l-(M p ))(l-(M h )) (60) 

= -2S W (e(h - 1) + f) + G^/(l-(M fc »(l-(M v » 
They can be seen as a simplification of the SCRPAx equations when the limit 

(QiQj) - 0, (Q l Q j ) - 0. (61) 

is assumed valid. At the same time, we see that the correlation energy can 
only be obtained using the expression for E^fr. 
Taking additionally the limit 

(Mi) -> 0, (62) 
the standard RPA matrices for the picket fence model 

Kv> A = 2S pp> (<P ~ l) + f ) - G » 

B ph A = G , (63) 

C% A = -26 hh , (e(h - \) + f ) + G , 

are recovered. 

7 Sum rules 

Given that the present RPA formalism is number conserving, we know that 
N\RPA) = Q\RPA), with N = ^ N h + ^ N p , (64) 

h P 

which satisfy 

n/2 n/2 n/2 

(AT) = 5^(2 - (M A » + £(M P > = n + £((M p=l > - (Mh=i)) = n (65) 
h=i p=i i=i 

due to (SB). 
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Given that 

Q 2 = (N 2 ) = ((Q + ]T(M p=i - M h=i ))(Q + £(M P= , - M h=j ))) 

i j 

= n 2 + 2ft ]T(M p=i - M h=i ) (66) 

i 

+ Y,(( M P^ M P=j) + {M h=i M h=j ) - (M p=i M h=j ) - (M h=i M p=j )) 

ij 

we get an additional condition over the expectation values, which is 

Y,(M p M pl )=Y,(M h M p ). (67) 

pp' hp 

This sum rule plays an important role in the discussions of the following 
section. 

8 The SCRPA ground state 

In all our previous derivations we always made use of the vacuum condition 
( |28|) for the ground state \SCRPA) of the present theory. However, at no step 
we have constructed the ground state explicitly. This situation is in fact quit 
common. Also Hartree-Fock and Hartree-Fock-Bogoliubov theories can be 
derived using analogous vacuum conditions, without constructing the ground 
state wave function explicitly (in spite of the fact that their construction is 
perfectly possible). The question then naturally arises whether from ( p8[ ) it 
is possible to explicitly construct the SCRPA ground state. 

To this purpose let us first consider the Q = 2 case. In this case there is 
only one particle state (p) and one hole state (h). There are two unperturbed 
(basis) states with 2 particles. They are 

\h) = \HF) and \p) = QlQ{\HF). (68) 

The first state has the two particles in the lower state (h), the second in the 
upper state (p). 

The SCRPA ground state can be expanded in this basis as 

| SCRPA) = a h \h)+ a p \p) with a\ + al = 1. (69) 
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The coefficients dh,a p must be determined from the vacuum condition (pffi), 
which in this simple case reads 

A \SCRPA) = (X p Ql - Y h Q h ) \SCRPA) = 0. (70) 

In this case X p = Xh = X and Y p = Yh = Y, and the removal operator R 
acting on the SCRPA vacuum yields the same equation, which is 

Xa p -Ya h = . (71) 

Together with the normalization condition, they have the solutions 

X Y 
VX 2 + Y 2 ' ap ~ VX 2 + Y 2 



a h = 7= : = , a p = . (72) 



In this way the SCRPA ground state is entirely determined by the RPA 
amplitudes X,Y. In the following section it will in fact be shown that the 
SCRPA can reproduces exactly the correlation and excitation energies in the 
case of Q — 2. 

Let us now consider the Q = 4 case. The phonon operators for the 
addition mode are 

4 = E Ql-Y. y h Q h = il x;q\ - ]T yjpQ fc> (73) 

p—l h=l p=l h=l 

were we have introduced the shorthand notation 



X£ = Xg/yJl-(M p ), Y£ = Y£/y/l-(M h ). (74) 
The SCRPA vacuum for 4 particles states has the form 
\SCRPA) = a \HF) + Y / ^ p Q{Q ] p \HF) + a 2 Q{ 1 Q{ 2 Ql 1 Ql 2 \HF). (75) 

hp 

The six coefficients do, a>hpi a 2 must be determined using the normalization 
and the vacuum conditions. They yield to the equations 

£^hp-n M ao = , Xga 2 -J2Y£a hp = . (76) 
p h 

There are four of these equations (2 for h = 1, 2 and two for p = 1,2) for 
each of the 2 possible /x's, giving a total of 8 equations, plus the normalization 
condition, to obtain the 6 a's. 
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These equations are consistent with the simplest HF (G=0) limit. In this 



case 



X£ = <5 W , Y£ = 0, (i,p,h = 1,2 (77) 

From Eq. ( j76|) we get ah p = for all h,p, and ci2 = 0. The only remaining 
variable to be determined is eto = 1.0 from the normalization condition. The 
HF ground state is recovered as expected. 

Given that the X£ 's are always finite, we can obtain four different ex- 
pressions for 02- They are 

Y 2 l Yl K, 1 

fl 2 - ^T a n + ^T fl 2i - 7-^12 + "-^22 
A 2 A 2 Aj 



2 



^ ^l 2 

^ fl U + ^i fl 21 = ^i«12 + ^i«22 (7 
A 2 A 2 aj Aj 



They imply 













Xj 2 ' 
















Af ' 


x* 


Xf 



(79) 

In general there is no reason for relations ( |79|) to be fulfilled. Indeed, in 
numerical examples they are not satisfied for any finite value of G. For larger 
spaces the situation does not improve. It implies that besides for Q = 2, the 
SCRPA ground state defined by Eq. fl28|) does not in fact exist. However, the 
theory is perfectly consistent in itself and ( p8|) must therefore be interpreted 
as an auxiliary relation which allows to express expectation values such as in 
(|29|) and ( |30|) in a well defined way in terms of the RPA amplitudes. 

In fact, using for example the Green's function method, as described 

in 



|13|, one can obtain expressions (|9],[3(]) via the residua of the Green's 
functions without ever using the vacuum condition ( p8| ) explicitly. On the 
other hand, the non-existence of \SCRPA) entails that we are not assured 
that our theory yields an upper bound of the ground state energy as in a truly 
Raleigh-Ritz variational principle. It also implies that our theory eventually 
contains some violation of the Pauli principle, in spite of the fact that we 
never made use of any boson approximation and always fully respected the 
commutation laws of the Fermi pair operators. In the next section, where 
we present the numerical results, we will return and further discuss some 
consequences of the non-existence of the SCRPA vacuum. 
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9 Numerical solutions 



In the present section the picket fence model is solved using the fully Self- 
Consistent RPA for different number of levels f2. The ground state and 
excited states energies, as well as the particle numbers and correlation func- 
tions, are compared with the exact ones, and with those obtained from the 
SCRPAi, the r-RPA and the RPA. For Vt = 2 analytical results are presented. 
In all the cases we fix the energy scale by setting e = 1. 

9.1 The n = 2 case 

This case was introduced in the previous section.. The Fermi level is at the 
energy 

A-f-f (80) 

Exact correlated states are built as linear combinations of the states basis 
states \h) and \p): 

\a) = c h \h) + Cp\p) with c 2 p + c 2 h = l, (81) 
which satisfy the eigenvalue equation 



with solutions 



E± = ±Vl + G 2 , c h = 



G E ± + l 



y/cp + (E ± + iy' ° p ] /G 2 + (E ± + l) 2 ' 

(83) 

The exact ground state |0) has energy E_ = - VI + G 2 . The correlation 
energy E COTT is 



^corr = (0|#|0) - (HF\H\HF) = -y/l + G 2 + 1 . (84) 
The exact particle number expectation values are 

(0|A4|0) = (0|M P |0) =2o 2 , (85) 
and the exact correlation functions are 

(0\M h M h \0) = (0\M p M p \0) = (0\M p M h \0) = Aaj (86) 
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The addition mode has 4 particles and only one state 

|4> = Q\\HF) , (87) 

with 

(4|#|4) = . (88) 
The exact excitation energy associated to this mode is 

E exc = (4\H\A) - (0\H\0) = = VTTCP. (89) 

The SCRPA equations for the addition mode are 



1-2GXY G(1 + 2F 2 ) I _ X\ 
-G(1 + 2Y 2 ) -1 + 2GXY [ Y ) ~ [ Y J iyUj 



with solutions 



G „ 1-E 1 



E 1>2 = ±VTT&, X= Y 

y/GP-(l-E 1>2 )* ^G 2 - (1 - £ lj2 ) 2 

(91) 

The addition mode has excitation energy E ± , which is exactly the same as 
-^exc deduced above, although the expressions for X, Y are quite different 
from those obtained for ai, a 2 - 
The occupation numbers are 

which again reproduce the exact results. 

The standard RPA equations for the addition mode are 

( —G -1 ) ( Y ) - *™ ( Y ) < 93 » 

with solutions 

E RPA = ±VT r G 2 (94) 

The RPA excitation energy exhibits the well known "collapse" when G — 1, 
and is a decreasing function of G, while the exact energy is an increasing 
function. 

The renormalized RPA (r-RPA) equations for the addition mode are 
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G 


exact 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.2 


-0.0198 


-0.0202 


-0.0198 


-0.0198 


-0.0198 


0.4 


-0.0770 


-0.0835 


-0.0770 


-0.0755 


-0.0770 


0.6 


-0.1661 


-0.2000 


-0.1362 


-0.1568 


-0.1661 


0.8 


-0.2806 


-0.4000 


-0.2065 


-0.2532 


-0.2806 


0.90 


-0.3454 


-0.5641 


-0.2415 


-0.3050 


-0.3454 



Table 1: Correlation energies as function of G for SI = 2 obtained with an 
exact calculation, with the RPA, r-RPA, SCRPAi and SCRPA methods. 



G 


exact 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.0 


1.0 


1.0 


1.0 


1.0 


1.0 


0.2 


1.0198 


0.9798 


0.9845 


1.0045 


1.0198 


0.4 


1.0770 


0.9165 


0.9572 


1.0139 


1.0770 


0.6 


1.1662 


0.8000 


0.9397 


1.0227 


1.1662 


0.8 


1.2806 


0.6000 


0.9379 


1.4155 


1.2806 


0.90 


1.3454 


0.4359 


0.9419 


1.5013 


1.3454 



Table 2: Energies of the first excited state as function of G for Vt = 2 ob- 
tained with an exact calculation, with the RPA, r-RPA, SCRPAi and SCRPA 
methods. 



1 + G 



2Y 2 



1+2Y 2 





E r - RPA [ v . (95) 



I n i _i _ 

\ ^X+2Y 2 L ^ 1+2Y 2 

Its solutions are not analytically simple. Numerical results can be seen in 
Tables [I] and ^|, together with the energies obtained with the other formalisms 
and the exact solutions. 

The contrast with the standard RPA exhibits the power of the fully self 
consistent RPA formalism, which in the Q = 2 case was shown to reproduce 
the exact results. 
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9.2 The Q = 4 case 

We label the four levels, from bottom to top, as h 2 , hi,px,p 2 . The exact 
eigenstates for 4 particles states, including the ground state, are 



a 



) = $\HF) + Y,<&QlQl\m + <5QlQlQlQl\HF)- (96) 

hp 



The coefficients cf. fulfill the normalization condition 

Ei<?r =i - ( 97 ) 

and are obtained as the eigen vectors of 



/ 


-4 


-G 


-G 


-G 


-G 





\ 


/ 


r a 
6 


\ 




( 


,,a 
< 


\ 




—G 


-2 


-G 


-G 





-G 






r a 








u h ipi 






-G 


-G 








—G 


-G 










= E a 




r a 

u h lP2 






-G 


-G 








-G 


-G 






r a 








r a 

U h 2 pi 






-G 





-G 


-G 


2 


-G 






r a 

^h 2 p 2 








r a 

^h 2 p2 




V 





-G 


-G 


-G 


-G 


4 


/ 


V 


r a 
c 2 


J 




V 


1 2 


J 



(98) 



The particle-hole symmetry implies that there are two degenerate unper- 
turbed states, which are 

QlQlJHF) and Q{ 2 Ql\HF), (99) 

both having zero unperturbed energies. They always contribute on equal 
footing to any correlated state, i.e. 

ctp2 = c a h2Pl . (100) 

The exact occupation numbers are 

(a\M h \a) = 2 (\c a hp f + \c a hp f + \c a 2 \ 2 ) , 

(a\M p \a) = 2 (\c a hip \ 2 + \c a h J 2 + \c a 2 \ 2 ) . (101) 

which due to ( |100[ ) imply (a\M h \a) = (a\M p=h \a); and 

(a\M h M h \a) = 2 (a|M ft |a), (a\M p M p \a) = 2 (a\M p \a), 

(a\M hl M h2 \a) = 4|c^| 2 , (a\M Pl M P2 \a) = 4 |cf | 2 , (102) 
(a\M p M h \a)=A\c a hp \ 2 + A\c^\ 2 . 
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G 


exact 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


-0.0123 


-0.0124 


-0.0123 


-0.0124 


-0.0124 


0.2 


-0.0520 


-0.0538 


-0.0510 


-0.0523 


-0.0524 


0.3 


-0.1225 


-0.1343 


-0.1158 


-0.1233 


-0.1248 


0.4 


-0.2264 


-0.2764 


-0.2009 


-0.2259 


-0.2334 


0.5 


-0.3644 


-0.5662 


-0.2961 


-0.3575 


-0.3795 



Table 3: Correlation energies as function of G for Q = 4 obtained with exact 
calculation, with the RPA, r-RPA, SCRPAi and SCRPA methods. 



The exact eigenstates for states with 6 particles (addition mode) are 



J24qI\hf) + J24q[q1 1 q1 2 \hf). 



(103) 



The coefficients c£ h are obtained as the eigen vectors of 



-3 


-G 


G 


G 


—G 


-1 


G 


G 


G 


G 


1 


-G 


G 


G 


—G 


3 



/ 



f 



E h 



(104) 



The exact excitation Eq XC energies for the addition mode are measured 
against the ground state energy, i.e. 



^exc 



Er — E a= 



(105) 



We have calculated the exact solutions, and compared them with the so- 
lutions of the SCRPA equations. This involved to solve simultaneously the 
SCRPA eigenvalue problem, Eqs. ([T9|) and to evaluate the number correla- 
tions by solving Eq. ([11]). In the following tables we present the correlation 
energies and the excitation energies of the first and second addition modes, 
obtained using four different methods: the exact results, the RPA ones, the 
SCRPA using (M^) = (M i )(M i ) (denoted SCRPA^ and the fully Self- 
Consistent RPA (SCRPA) which refers to the treatment described above. 

In the three tables ^],£|,|5] the improvement offered by the SCRPA as com- 
pared with the RPA is clearly seen. Both for the correlation and the excita- 
tion energies the RPA gives only a crude approximation, while the SCRPA 
results are very close to the exact ones. 
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G exact RPA r-RPA SCRPA! SCRPA 



0.0 


1.0 


1.0 


1.0 


1.0 


1.0 


0.1 


1.0027 


0.9864 


0.9872 


1.0041 


1.0027 


0.2 


1.0115 


0.9402 


0.9494 


1.0233 


1.0126 


0.3 


1.0276 


0.8493 


0.8917 


1.0648 


1.0348 


0.4 


1.0520 


0.6898 


0.8248 


1.1279 


1.0734 


0.5 


1.0853 


0.3745 


0.7603 


1.2065 


1.1274 



Table 4: Excitation energy of the first addition mode as function of G for 
f2 = 4 obtained with exact calculation, with the RPA, r-RPA, SCRPAi and 
SCRPA methods. 



G 


exact 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.0 


3.0 


3.0 


3.0 


3.0 


3.0 


0.1 


3.0104 


3.0012 


3.0014 


3.0108 


3.0104 


0.2 


3.0465 


3.0060 


3.0075 


3.0499 


3.0468 


0.3 


3.1156 


3.0164 


3.0212 


3.1259 


3.1168 


0.4 


3.2237 


3.0338 


3.0446 


3.2422 


3.2256 


0.5 


3.3738 


3.0594 


3.0784 


3.3967 


3.3784 



Table 5: Excitation energy of the second addition mode as function of G for 
Q = 4 obtained with exact calculation, with the RPA, r-RPA, SCRPAx and 
SCRPA methods. 
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0.12 



exact — 



0.1 - 




0.08 - 



>/2 



0.06 - 



0.04 - 



0.02 - 











0.1 



0.2 



0.3 



0.4 



0.5 



G 



Figure 1: Occupation numbers (M p )/2 as a function of G for Q — 4. 

The correlation energies listed in Table are slightly better described 
with the SCRPAx than with the SCRPA. It is also noticeably that both 
the RPA and SCRPA correlation energies are lower than the exact ones, 
but the SCRPA! energies are slightly higher. The situation is very different 
when excitation energies are analyzed (Tables [| and ^). While again the 
RPA energies are lower than the exact ones, both the SCRPA and SCRPAi 
energies are higher. As mentioned above, they are very close to the exact 
ones, but in this case the SCRPA ones are clearly the best, reproducing with 
high accuracy the exact results. 

To obtain a deeper understanding of the above discussed results a study 
of the occupation numbers in the different approaches is in order. They are 
plotted in Fig. 1 as a function of G. The occupation numbers (M p )/2 obtained 
using the exact solutions are represented by thin lines, those obtained using 
the RPA by small dots, the SCRPAi by thick lines and the SCRPA by large 
dots. Upper curves display the occupations for p=l, lower curves for p=2. 

As it is well known, RPA calculations predict far more correlations than 
they are found in the exact ground state. For both states this is evident 
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a j 




Exact solutions 




i 


j 


(M h M hj ) 


(M hl M Pj ) 


\ M Pi M Pi) 


1 


1 


0.2556 


0.1967 


U.zooo 


i 
i 


o 
Z 


0.0028 


0.0616 


U.UUZo 


2 


1 


0.0028 


0.0616 


0.0028 


2 


2 


0.0864 


0.0275 


0.0864 


bj 




SCRPA solutions 




i 


j 


(M hi M hi ) 


(M hi M Pi ) 


(M Pi M Pi ) 


1 


1 


0.2251 


0.1774 


0.2251 


1 


2 


0.0033 


0.0603 


0.0033 


2 


1 


0.0033 


0.0603 


0.0033 


2 


2 


0.0816 


0.0281 


0.0816 


c) 




product matrix 


i 


j 


(M hi ){M hi ) 


(M hi ){M Pi ) (M Pt )(M P] 


1 


1 


0.0127 


0.0127 


0.0127 


1 


2 


0.0046 


0.0046 


0.0046 


2 


1 


0.0046 


0.0046 


0.0046 


2 


2 


0.0017 


0.0017 


0.0017 



Table 6: Correlation functions (MjMj) for G=0.4, Q — 4 . Insert a) shows 
the exact results, insert b) the SCRPA results and insert c) the product of 
single occupations 

here: the small dotted curves increase as function of G faster then the exact 
ones, and diverge close to the point of collapse of the RPA. On the other 
hand, we learn from Fig. 1 that the SCRPA under predict the ground state 
correlations. The SCRPAi occupations deviates from the exact ones nearly 
as much as the RPA occupations, but underestimating them. It is remark- 
able that the SCRPAi energies are so good, given these discrepancies. The 
SCRPA occupations are very close to the exact ones. 

The correlation functions (MjMA are shown in Table |]for G=0.4. 

As can be seen in Table |5] the occupation number correlations obtained 
self-consistently from the SCRPA equations are very close to the exact ones, 
even for a value of G as large as 0.4. On the other hand the approximation 
(MiMj) pa (Mi)(Mj) underestimates the correlations. It could explain the 
small occupation numbers obtained when this approximation is used in the 
SCRPA equations. 
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E(M P M P ,) J2(M h M p ) 

vvL t}E 

exact 0.3476 0.3476 

SCRPA 0.3133 0.3260 
Table 7: Sum rule fl67D for the exact and SCRPA calculations for G=0.4, 

n = 4 

Although the SCRPA results presented here are in general extremely close 
to the exact results, considering correlation or excitation energies, occupation 
or correlation numbers, there is a caveat which must be mentioned here. It 
concerns the sum rule (|5T[) . 

The sums of correlation numbers obtained using the exact calculation and 
the SCRPA one are very similar. They are shown in Table [7|. However, while 
the exact results fulfill the sum rule fl6T|), the SCRPA do not fulfill this sum 
rule. This shortcoming certainly stems from the fact that, as we explained 
in Section 8, the SCRPA ground state as defined in Eq. fl28|) only exists 
for the special case Q = 2. In all the other cases the use of relation (p8|) 
implies some degree of approximation entailing, for example, the violation of 
the above sum rule. 



9.3 The Q = 10 case 

The results for Q = 10 repeat the general patterns found for Q = 4. In Table 
| the ground state correlation energies are presented as a function of the 



pairing strength G, in Tables and 10 the excitation energies of the first and 



second addition mode in a system with N=12 fermions, measured relative to 
the ground state of the system with N=10 are shown, respectively. 

The SCRPA does not show any collapse, and the energies are far closer 
to the exact energies than the RPA ones. For the ground state correlation 
energy the ones obtained with the SCRPAi formalism are slightly closer to 
the exact ones than the SCRPA energies, while the opposite is true for the 
first and second excited states. For these excitation energies the SCRPA 
reproduces with very high accuracy the exact values. It is worth to mention 
that the exact excitation energies of the first and second addition modes 
monotonically increase as function of G. This behavior is well reproduced 
by the SCRPA, while in the standard RPA these excitation energies are 
decreasing. This implies that the screening of the bare interaction is so strong 
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G 


exact 


RPA 


SCRPAi 


SCRPA 


.00 


0.0000 


0.0000 


0.0000 


0.0000 


.05 


-0.0086 


-0.0086 


-0.0086 


-0.0086 


.10 


-0.0364 


-0.0367 


-0.0365 


-0.0365 


.20 


-0.1669 


-0.1756 


-0.1686 


-0.1691 


.30 


-0.4350 


-0.5419 


-0.4424 


-0.4497 


.33 


-0.5505 


-0.8181 


-0.5594 


-0.5725 


.34 


-0.5931 


#** 


-0.6023 


-0.6180 


.35 


-0.6379 


*** 


-0.6473 


-0.6658 


.36 


-0.6850 


#*# 


-0.6943 


-0.7160 



Table 8: Correlation energies as function of G obtained with exact calcula- 
tion, with the RPA, SCRPAi and SCRPA methods, for ft = 10. 



that even the sign of the interaction is turned around, bringing the SCRPA 
solution closer to the exact one. This is a quite remarkable achievement of 
SCRPA. 

The sum rule (|67j) is again violated in the ft = 10 shown in Table 



11. It must be said that the difference between both sums, of the order 0.02, 



must be compared with N 2 = 100. The violation of the sum rule is therefore 
very small. On the other side, it exhibits the limitations of the best possible 
SCRPA treatment of this problem. 

9.4 The n = 24 and Q = 100 cases 

The correlation energies for Q = 24 are shown in Table ITI3. The excitation 



energies corresponding to the first and second addition modes are presented 
in Tables [H| and [14] respectively. All the energies are shown as functions 
of the pairing strength G given in the first columns. The second column 



exhibits the exact results in Table 12. RPA, r-RPA, SCRPAi and SCRPA 



energies are presented in the following four columns. 

The large Q limit is being approached, where all the RPA descriptions 
show their best performances. For the correlation energies even the standard 
RPA predictions are good, except for G = 0.25, close to the point of collapse. 
For the excitation energies both the RPA and r-RPA results decrease when 
G increases, while the SCRPA energies are increasing, as in previous cases, 
and in agreement with the exact energies (not shown). At the same time 
the differences between the SCRPAi and SCRPA energies are vanishing. It 
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G 


exact 


RPA 


SCRPAi 


SCRPA 


.00 


1.0000 


1.0000 


1.0000 


1.0000 


.05 


1.0003 


0.9940 


1.0005 


1.0003 


.10 


1.0011 


0.9732 


1.0034 


1.0014 


.20 


1.0053 


0.8604 


1.0279 


1.0119 


.30 


1.0143 


0.5257 


1.0970 


1.0539 


.33 


1.0184 


0.2574 


1.1266 


1.0758 


.34 


1.0199 


*** 


1.1372 


1.0840 


.35 


1.0216 


*** 


1.1481 


1.0927 


.36 


1.0233 


*** 


1.1592 


1.1018 



Table 9: Excitation energy of the first addition mode as function of G ob- 
tained with exact calculation, with the RPA, SCRPAi and SCRPA methods, 
for tt = 10. 



G 


exact 


RPA 


SCRPAi 


SCRPA 


.00 


3.0000 


3.0000 


3.0000 


3.0000 


.05 


3.0010 


2.9971 


3.0011 


3.0010 


.10 


3.0056 


2.9881 


3.0063 


3.0057 


.20 


3.0390 


2.9515 


3.0460 


3.0406 


.30 


3.1442 


2.8935 


3.1657 


3.1507 


.33 


3.1998 


2.8736 


3.2244 


3.2070 


.34 


3.2214 


*** 


3.2467 


3.2287 


.35 


3.2449 


*** 


3.2703 


3.2518 


.36 


3.2701 


*** 


3.2952 


3.2765 



Table 10: Excitation energy of the second addition mode as function of 
G obtained with exact calculation, with the RPA, SCRPAi and SCRPA 
methods, for Q = 10. 



J2(M p M pl ) J2(M h M p ) 

pp' hp 



SCRPA 0.4509 



0.4736 



Table 11: Sum rule (|6"7j ) for the exact and SCRPA calculations for f2 = 10 
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G 


exact 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.00 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.05 


-0.0218 


-0.0218 


-0.0218 


-0.0218 


-0.0218 


0.10 


-0.0953 


-0.0961 


-0.0955 


-0.0955 


-0.0955 


0.15 


-0.2373 


-0.2432 


-0.2358 


-0.2386 


-0.2388 


0.20 


-0.4736 


-0.5096 


-0.4623 


-0.4794 


-0.4809 


0.25 


-0.8452 


-1.1420 


-0.7808 


-0.8607 


-0.8694 



Table 12: Correlation energies as function of G obtained with exact calcula- 
tion, with the RPA, r-RPA, SCRPAi and SCRPA methods for Q = 24. 



strongly supports the use of the simpler SCRPAx approach for larger fi, 
where the numerical effort needed to solve the SCRPA equations in the fully 
self-consistent formalism is not justified. 

This is the case for Q = 100, whose energies are shown in Tables [15], [T^ 
and [17| as functions of the pairing strength G. Table [15] exhibits the exact 
correlation energies in the second columns, and the RPA, r-RPA and SCRPAx 
in the third, fourth and fifth columns, respectively. The agreement between 
the exact and SCRPA is excellent, even for G = 0.19, where the RPA has 
collapsed. 

Tables [16] and [L7| present the RPA, r-RPA and SCRPAi excitation en- 
ergies of the first and second addition modes as function of G. Again only 
the SCRPA energies are increasing when G increases, as is expected from 
the exact results. The latter are not shown, since it would have needed extra 
numerical effort without bringing new insight to the problem. It should be 
mentioned that Q = 100 constitutes an extremely large configuration space 
which by no means can be handled with ordinary diagonalization procedures. 

It is well known that the RPA collapse reflects a qualitative change in the 
mean field, which after the collapse is dominated by superfluid correlations 



24f . It must be remembered that, while the energetics and the occupation 
numbers obtained with the SCRPA are very close to the exact ones, the wave 
functions around and beyond the value of G at which standard RPA collapses, 
being far better than those obtained with RPA or r-RPA, can nonetheless 
have an overlap with the exact wave function of less than 50% p5] . In this 



case the SCRPA must be extended to the deformed, i.e. superfluid, basis 

inn. 
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G 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.00 


1.0000 


1.0000 


1.0000 


1.0000 


0.05 


0.9911 


0.9912 


1.0004 


1.0001 


0.10 


0.9574 


0.9594 


1.0035 


1.0011 


0.15 


0.8806 


0.8936 


1.0153 


1.0060 


0.20 


0.7155 


0.7818 


1.0458 


1.0233 


0.25 


0.2132 


0.6279 


1.1040 


1.0661 



Table 13: Excitation energy of the first addition mode as function of G 
obtained with the RPA, r-RPA, SCRPAi and SCRPA methods for = 24. 



G 


RPA 


r-RPA 


SCRPAi 


SCRPA 


0.00 


3.0000 


3.0000 


3.0000 


3.0000 


0.05 


2.9944 


2.9944 


3.0007 


3.0006 


0.10 


2.9755 


2.9761 


3.0057 


3.0049 


0.15 


2.9398 


2.9425 


3.0233 


3.0200 


0.20 


2.8836 


2.8926 


3.0705 


3.0622 


0.25 


2.8044 


2.8297 


3.1758 


3.1618 



Table 14: Excitation energy of the second addition mode as function of G 
obtained with the RPA, r-RPA, SCRPAi and SCRPA methods for Q = 24. 



G 


exact 


RPA 


r-RPA 


SCRPA 


0.00 


0.0000 


0.0000 


0.0000 


0.0000 


0.05 


-0.0946 


-0.0946 


-0.0946 


-0.0946 


0.10 


-0.4250 


-0.4274 


-0.4247 


-0.4256 


0.12 


-0.6467 


-0.6540 


-0.6454 


-0.6483 


0.14 


-0.9365 


-0.9575 


-0.9314 


-0.9405 


0.16 


-1.3127 


-1.3768 


-1.2942 


-1.3217 


0.17 


-1.5413 


-1.6628 


-1.5067 


-1.5543 


0.18 


-1.8027 


-2.0914 


-1.7395 


-1.8209 


0.19 


-2.1028 


*** 


*** 


-2.1273 



Table 15: Correlation energies as function of G obtained with exact calcula- 
tion, with the RPA, r-RPA and SCRPAi methods for = 100. 
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G 


RPA 


r-RPA 


SCRPAx 


0.00 


1.0000 


1.0000 


1.0000 


0.05 


0.9857 


0.9859 


1.0004 


0.10 


0.9210 


0.9252 


1.0056 


0.12 


0.8657 


0.8774 


1.0126 


0.14 


0.7762 


0.8087 


1.0262 


0.16 


0.6186 


0.7144 


1.0504 


0.17 


0.4829 


0.6575 


1.0678 


0.18 


0.2181 


0.5958 


1.0895 


0.19 


*** 


*** 


1.1154 



Table 16: Excitation energy of the first addition mode as function of G 
obtained with the RPA, r-RPA and SCRPAx methods for SI = 100. 



G RPA r-RPA SCRPAx 

0.00 3.0000 3.0000 3.0000 

0.05 2.9896 2.9897 3.0006 

0.10 2.9489 2.9499 3.0082 

0.12 2.9193 2.9218 3.0182 

0.14 2.8788 2.8843 3.0381 

0.16 2.8245 2.8362 3.0761 

0.17 2.7912 2.8080 3.1059 

0.18 2.7530 2.7773 3.1457 

0.19 *** *** 3.1985 

Table 17: Excitation energy of the second addition mode as function of G 
obtained with the RPA, r-RPA and SCRPAi methods for ft = 100. 
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10 Conclusions 



In this work we presented for the first time a fully Self- Consistent RPA 
treatment of the picket fence model. This involved developing equations for 
some correlation functions which had to be solved simultaneously with the 
SCRPA eigenvalue problem. 

One of the remarkable results of the present work is that for Q = 2, i.e. 
for the 2 particles case, SCRPA provides the exact energies of the problem. 
This is not at all trivial because usually many body approaches deteriorate 
as the number of particles decreases. For higher number of levels SCRPA is, 
of course, not any longer exact but the results for the correlation energy and 
the low lying part of the spectrum are in excellent agreement with the exact 
ones. 

In an earlier work we had used the factorization ansatz (NiNj) ~ (JVi) (Nj) 
for the density-density correlation functions. As mentioned already, we did 
not make any of these approximations in the present work. Nevertheless we 
could verify a posteriori that the above approximation yields still excellent 
results as compared with the exact solutions. This is more noticeable when Q 
is increased. Actually we found that the above factorization approximation 
also works well for other models [I5|. Therefore one may conjecture that it 



is quite generally valid. It is also quite useful, since it reduces the numerical 
work considerably. For example with its use we solved the SCRPA equations 
for Q = 100 without problems. 

In spite of the excellent results we could not confirm previous conjectures, 
at least for the ansatz used in this work, that SCRPA yields an upper limit 
to the ground state energy. Our values are very accurate but slightly below 
the exact values. This finding is certainly linked to the fact that for the 
present case a Self- Consistent RPA ground state wave function does not 
exist except for Q = 2, as we showed in Section 8. This most likely entails 
a violation of the Pauli principle which, although small, in turn gives raise 
to overbinding. This also explains the slight violation of the sum rule (67). 
One should, however, mention that the picket fence model, because of the 
low degeneracy of the levels, tests the Pauli principle in a most sensitive way. 
This is, for example, born out in the fact that the low lying excited states, as 
a function of G, raise instead of getting lower as one could expect from the 
attractive nature of the interaction. It is another outstanding result of the 
present application of the SCRPA theory that its results follow very closely 
this trend of the excited states implying that the renormalization of the bare 
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interaction due to the self-consistency effect can even turn around the sign 
of the interaction. 

We also should mention that we have treated only the non-superfluid 
phase. It is well known from our previous studies that for interaction strengths 
beyond the collapse of standard RPA one must switch to the 'deformed' or 
superfluid basis, otherwise the results deteriorate or the iteration procedure 
does not converge. The superfluid domain of the present model shall be 
studied in a future work. 

All in all we can conclude that in this first full application of SCRPA to 
a large scale problem the expectation we had gained from earlier studies in 
the performance of SCRPA has been fully confirmed. Indeed the results did 
not fail in any major quantitative respect. It may therefore seem worthwhile 
to push the applications of SCRPA to more realistic Hamiltonians. 
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Appendix: Some useful relations 

R\, Ry 

[M p , R x ] 
M h ,4 



— V vAyA' 1 ~ m p I Y X Y X ' 1 ~ M fe 
l^p I p I p i-(m p ) ^ l^h^h^h i-(M h ) ' 

vp. yn' l-M p yMyM' kzMh 

/^ p s\ p s\ p i_( Mp ) Z^h 1 h 1 h l-(M h ) ' 

= - Ep X p Y p X l-{M v ) + Eft Y h X hl-(M h ) ' ( 106 ) 

= -21?{E W *f 4i + £ Al Y*R Xl } 

2>T{E M1 ^4i + Eai x^ju 



yAiyA' 

P2 P2 



yAyA2 
Pi Pi 

^ (1 - <M Pl » (1 - (M P2 » 



(RxR\ 2 RxiRy) - 
((1 - M P1 )(1 - M p2 )) 



yAyA2 
_ pi pi 



A fti A fti 



p ^ (1 - (M P1 » (1 - (M fcl » 



A /ii A fti 



yAiyA' 
_ Pi Pi 

fa (1 - (M P1 » (1 - (M fcl » 

x^A "V^A2 



vM yA' 

A ft2 A ft2 



(M h2 )) 



((1 - M P1 )(1 - M hl )) 

((1-M M )(1-M fcl )> 
((l-MfcJCl-Mfc,)) 



(107) 



(i?A [-Rai , -Ra 2 



2E 



pi 



y Ai y- A 2 y A y A' 

Pi Pi Pi Pi 

(1 " (M P1 )) 



2E 

hi 



v^i v^z yA yA' 

A fej A fe! A fc 1 A fe 1 

^ (l-(M fcl » 
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yAiyA 2 
I Pi Pi _ 

tp, (1 - < M P1» (1 - < M P 2 » 



yAyA' 

((1-M pi )(l-M p2 )> 



yAiyA 2 yA yA' 

yAyA' yAi t^A 2 

_ y pi pi fei fe \ x ((i-M Pl )(i -M ftl )) 
P V 1 (i-(^ 1 ))(i-(^ 1 » U PlA J/ 



yA yA' 



+ £(l-<M fcl »(l-(M fca » 



((1 - M hl )(l - M h2 )) 



(108) 



(R\R\ l R X2 R y ) — (R\ R\ 1 ,R\ 2 Ry) + (R\R Xl R\ 2 R x ,) 

yAiyA2yAyA' x^Ai -^A 2 t^A yA' 

_ 2 y^ pi pi pi pi _ ^/VVHi 

yA 2 yA' /yAyAi , yA x yA\ 

+ X! % P / ,/\vi u (C 1 " M pi)( X " M p 2 )) 



PlP2 



(1-(M P1 »(1-(M P2 » 



fyA 2 yA' i yA' vA2VyAi yA _i_ yA yAlN 

^ 1 Pl ^TVTT^ 77J u Pl — U 1 -^)^-^)) ( 109 ) 



P\h\ 



(l-(M pl »(l-(M fcl » 

yA 2 yA' / yAi x^A I yA yAi\ 



(1 - (M hl ))(l - (M h2 )) 



(A* An ^2^1') — 

yMiyM2yMy^' Y/u y^2 y/i \V 

£ (l-(M fcl » ~ tT (I-(^pi)) 

yM2y£t' /yM yMi _i_ yMiyMA 

-g, ^a-lg^S))^ } <(1 " MJ(1 - (110) 
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